Molecular dynamics simulation of isothermal crystallisation of polymer chains around single polymer lamella

The isothermal crystallisation of polyethylene (PE) chains around single PE lamella in vacuum is investigated by molecular dynamic simulation. The crystallisation process is analysed in terms of the orientational order parameters, principal moments of inertia for the simulated systems. The effects of charge interactions between the polymer chains and lamella are discussed. It is found that the crystallisation process for uncharged systems can be divided into three stages: (1) adsorption, (2) orientation and (3) arrangement. The single polymer lamella changes a little during the three stages. PE chains are arranged parallel to the chain direction of the stems in the crystalline state. When considering the effect of charge interactions between the polymer chains and lamella, a different crystallisation process appears. The single polymer lamella is affected by the charged polymer chains.     

How good are comparative models in the understanding of protein dynamics?

The 3D structure of a protein is essential to understand protein dynamics. If experimentally determined structure is unavailable, comparative models could be used to infer dynamics. However, the effectiveness of comparative models, compared to experimental structures, in inferring dynamics is not clear. To address this, we compared dynamics features of ~800 comparative models with their crystal structures using normal mode analysis. Average similarity in magnitude, direction, and correlation of residue motions is >0.8 (where value 1 is identical) indicating that the dynamics of models and crystal structures are highly similar. Accuracy of 3D structure and dynamics is significantly higher for models built on multiple and/or high sequence identity templates (>40%). Three-dimensional (3D) structure and residue fluctuations of models are closer to that of crystal structures than to templates (TM score 0.9 vs 0.7 and square inner product 0.92 vs 0.88). Furthermore, long-range molecular dynamics simulations on comparative models of RNase 1 and Angiogenin showed significant differences in the conformational sampling of conserved active-site residues that characterize differences in their activity levels. Similar analyses on two EGFR kinase variant models highlight the effect of mutations on the functional state-specific αC helix motions and these results corroborate with the previous experimental observations. Thus, our study adds confidence to the use of comparative models in understanding protein dynamics.     

Analysis of non-peptidic compounds as potential malarial inhibitors against Plasmodial cysteine proteases via integrated virtual screening workflow

Falcipain-2 (FP-2) and falcipain-3 (FP-3), haemoglobin-degrading enzymes in Plasmodium falciparum, are validated drug targets for the development of effective inhibitors against malaria. However, no commercial drug-targeting falcipains has been developed despite their central role in the life cycle of the parasites. In this work, in silico approaches are used to identify key structural elements that control the binding and selectivity of a diverse set of non-peptidic compounds onto FP-2, FP-3 and homologues from other Plasmodium species as well as human cathepsins. Hotspot residues and the underlying non-covalent interactions, important for the binding of ligands, are identified by interaction fingerprint analysis between the proteases and 2-cyanopyridine derivatives (best hits). It is observed that the size and chemical type of substituent groups within 2-cyanopyridine derivatives determine the strength of protein–ligand interactions. This research presents novel results that can further be exploited in the structure-based molecular-guided design of more potent antimalarial drugs.     

Hard-sphere/pseudo-particle modelling (HS-PPM) for efficient and scalable molecular simulation of dilute gaseous flow and transport

Continuum methods are not accurate enough for flows at high Knudsen numbers, whereas rigorous molecular dynamics (MD) methods are too costly for simulations at practical dimensions. Hard-sphere (HS) model is a simplified MD method efficient for dilute gaseous flow but is of poor parallelism due to its event-driven nature, which sets a strong limitation to its large-scale applications. In this work, pseudo-particle modelling, a time-driven modelling approach is coupled with HS model to construct a scalable parallel method capable of simulating flows and transport processes at high Knudsen numbers without losing necessary molecular details in describing their macro-scale behaviours. The method is validated in several classical simulation cases and its performance is evaluated to be favourable. To demonstrate the potential applications of this method, we also simulate the diffusion of small molecules in multi-scale porous media which is related to catalysis, material preparation and micro chemical engineering in the long term.     

Long-term changes in a Wisconsin Fagus-Acer forest in relation to glaze storm disturbance

Changes in the composition of a Fagus-Acer (Beech-Sugar maple) forest in southeastern Wisconsin over a 16-yr period from 1971 to 1987 are analyzed in relation to a severe glaze (ice) storm disturbance occurring within the census period. Landscape topography created ‘windward’ and ‘leeward’ forest aspects with respect to storm severity, which resulted in greater canopy opening on the windward aspect. In the tree stratum, most species remained stable in density and most of the common species increased in basal area into larger size classes. However, Fagus grandifolia, Ulmus rubra, and the small tree Ostrya virginiana suffered net losses that suggest synergistic effects between glaze storm disturbance and other factors upon tree mortality. In the sapling stratum, canopy opening strongly promoted release of shade-to levant Acer sac-charum. On the windward forest aspect, sapling densities of less shade-tolerant species also increased, in contrast to the absence of such increases on the leeward forest aspect. In the shrub (regeneration) stratum, species responses were heterogeneous. Regeneration of most species increased over the 16-yr period, and some less shade-tolerant species showed increased regeneration differentially on the windward forest aspect. Overall, disturbance appears to have accelerated forest succession toward increased dominance by A. saccharum and persistence of both Fagus and Tilia americana through their capacities for root sprouting. However, forest succession was retarded somewhat on the windward aspect through increased recruitment of less-shade tolerant species. These results parallel those of other studies of glaze storm disturbance, and they illustrate how spatially heterogeneous disturbance intensity may contribute to maintenance of forest diversity.     

Atomic insight into prion disorder: An intricate detail gained by 0.5 μs molecular dynamics simulation of preventive G127V and deleterious D178V mutation in prion protein

In this study we are looking into two contradicting mutations found in prion protein (PrP) viz G127V and D178V, that are reportedly protective and pathogenic, respectively. Despite significant advances in comprehension of the role of pathogenic mutations, the role of protective mutation in amyloid fold inhibition still lacks a substantial basis. To understand the structural basis of protective mutation, molecular dynamics simulation coupled with protein-protein docking and molecular mechanics/Poisson-Boltzmann surface area analysis was used to understand the instant structural variability brought about by these mutations alone and in combination on PrP and prion-prion complex. Atomic-scale investigations successfully revealed that the binding pattern of prion-prion varies differentially in protective and pathogenic mutations with secondary structure showing distinct contrasting patterns, which could supposedly be a critical factor for differential prion behavior in protective and pathogenic mutations. Considering the reported role of an amyloid fold in prion-prion binding, the contrasting pattern has given us a lead in comprehending the role of these mutations and has been used in this study to look for small molecules that can inhibit amyloid fold for prion-prion interaction in pathogenic mutant carrying PrP.     

Identification of the phospholipid binding regions of the envelope E protein of flaviviruses by molecular dynamics

Abstract

The Flavivirus genus comprise several important human pathogens, including dengue, West Nile, Yellow fever, Japanese encephalitis, Zika, and tick-borne encephalitis viruses. These enveloped viruses affect more than 2 billion people in the world, mainly in less developed countries. Although some vaccines exist for some flaviviruses, these vaccines are not universally available due to many factors and since their infections are a world-wide public health issue, the development of antiviral molecules is fundamental. Flavivirus membranes, through the help of the envelope E glycoprotein, fuse with endosomal compartments in a pH-dependent way to release their genome into the cytoplasm and require specific lipids, such as bis(monoacylglycero)phosphate (BMP), for efficient fusion. The fundamental role the envelope E protein has on viral entry and membrane fusion suggest that it is an essential antiviral target. In this work, we have used atomistic molecular dynamics simulations to study the binding of the head-group of BMP to the tip of the envelope E proteins of ZIKV, DENV, TBEV and JEV viruses whose three-dimensional structures are known. Our results indicate that, apart from the fusion loop, there are different amino acid residues in different regions of the envelope E proteins of flaviviruses capable of binding the head-group of BMP. These regions should work together to accomplish the binding and fusion of the envelope and endosomal membranes and represent a new target to develop and design potent and effective antiviral agents capable of blocking flavivirus-endosome membrane fusion.      

Interactions of a Bacterial RND Transporter with a Transmembrane Small Protein in a Lipid Environment

1. Download : Download high-res image (350KB)
2. Download : Download full-size image

     

Molecular docking and simulation of Curcumin with Geranylgeranyl Transferase1 (GGTase1) and Farnesyl Transferase (FTase)

Protein prenylation is a posttranslational modification that is indispensable for translocation of membrane GTPases like Ras, Rho, Ras etc. Proteins of Ras family undergo farnesylation by FTase while Rho family goes through geranylgeranylation by GGTase1. There is only an infinitesimal difference in signal recognition between FTase and GGTase1. FTase inhibitors mostly end up selecting the cells with mutated Ras proteins that have acquired affinity towards GGTase1 in cancer microcosms. Therefore, it is of interest to identify GGTase1 and FTase dual inhibitors using the docking tool AutoDock Vina. Docking data show that curcumin (from turmeric) has higher binding affinity to GGTase1 than that of established peptidomimetic GGTase1 inhibitors (GGTI) such as GGTI-297, GGTI-298, CHEMBL525185. Curcumin also interacts with FTase with binding energy comparable to co-crystalized compound 2-[3-(3-ethyl-1-methyl-2-oxo-azepan-3-yl)-phenoxy]-4-[1-amino-1-(1-methyl-1h-imidizol-5-yl)-ethyl]-benzonitrile (BNE). The docked complex was further simulated for 10 ns using molecular dynamics simulation for stability. Thus, the molecular basis for curcumin binding to GGTase1 and FTase is reported.     

De novo inference of protein function from coarse‐grained dynamics

Inference of molecular function of proteins is the fundamental task in the quest for understanding cellular processes. The task is getting increasingly difficult with thousands of new proteins discovered each day. The difficulty arises primarily due to lack of high‐throughput experimental technique for assessing protein molecular function, a lacunae that computational approaches are trying hard to fill. The latter too faces a major bottleneck in absence of clear evidence based on evolutionary information. Here we propose a de novo approach to annotate protein molecular function through structural dynamics match for a pair of segments from two dissimilar proteins, which may share even <10% sequence identity. To screen these matches, corresponding 1 µs coarse‐grained (CG) molecular dynamics trajectories were used to compute normalized root‐mean‐square‐fluctuation graphs and select mobile segments, which were, thereafter, matched for all pairs using unweighted three‐dimensional autocorrelation vectors. Our in‐house custom‐built forcefield (FF), extensively validated against dynamics information obtained from experimental nuclear magnetic resonance data, was specifically used to generate the CG dynamics trajectories. The test for correspondence of dynamics‐signature of protein segments and function revealed 87% true positive rate and 93.5% true negative rate, on a dataset of 60 experimentally validated proteins, including moonlighting proteins and those with novel functional motifs. A random test against 315 unique fold/function proteins for a negative test gave >99% true recall. A blind prediction on a novel protein appears consistent with additional evidences retrieved therein. This is the first proof‐of‐principle of generalized use of structural dynamics for inferring protein molecular function leveraging our custom‐made CG FF, useful to all. Proteins 2014; 82:2443–2454. © 2014 Wiley Periodicals, Inc.